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Kurthia massiliensis strain JC30^ sp. nov. is the type strain of K. massiliensis sp. nov., a new 
species within the genus Kurthia. This strain, whose genome is described here, was isolated 
from the fecal flora of a healthy patient. K. massiliensis is a Gram-positive aerobic rod. Here 
we describe the features of this organism, together with the complete genome sequence and 
annotation. The 3,199,090 bp long genome contains 3,240 protein-coding genes and 86 
RNA genes, including between 3 and 4 rRNA genes. 



Abbreviations: EMBL: European Molecular Biology Laboratory 



Introduction 

Kurthia massiliensis strain ]C30^ [CSUR 141^ = 
DSM 24639^) is the type strain of K. massiliensis 
sp. nov. This bacterium is a Gram-positive, strictly 
aerobic rod that is capsulated, and motile by 
peritrichous flagella. This organism w^as originally 
isolated from the stool of a healthy Senegalese pa- 
tient as part of a "culturomics" study aimed at cul- 
tivating all species w^ithin human feces, individual- 

ly- 

Currently, "the gold standard" for defining bacte- 
rial species is DNA-DNA hybridization [1]. But this 
method is time-consuming and the inter- 
laboratory reproducibility is poor. Fortunately, 
the development of PGR and next-generation se- 
quencing technologies have led to reliable and re- 
producible 16S rRNA comparison methods w^ith 
generally agreed upon cutoff values that enable 
the taxonomic classification of new^ species for 
many bacterial genera [2]. To describe new^ bacte- 
rial taxa, the use of a polyphasic approach v\^as 
proposed [3] that includes their genome sequence, 
MALDI-TOF spectrum and main phenotypic char- 
acteristics [habitat. Gram-stain reaction, cultiva- 
tion conditions, cell w^all structure and metabolic 
characteristics). 

The genus Kurthia w^as created in 1885 by 
Trevisan [4] in honor of Kurth w^ho described the 
first species. Bacterium zopfii, isolated from the 
intestinal contents of chickens. As the stool sam- 
ples had been stored at room temperature and the 



bacteria v\^ere strictly aerobic, it v\^as assumed that 
the samples w^ere contaminated by Kurthia, w^hich 
multiplied during storage. The name Kurthia was 
first published in the seventh edition of Bergey's 
Manual of Determinative Bacteriology [5] and w^as 
included in the Approved Lists of Bacterial Names 
[6]. Currently, Kurthia includes 3 species: K zopfii, 
K gibsonii [7] and K sibirica [8]. The bacteria are 
members of the phylum Firmicutes, and the family 
Planococcaceae. There is no evidence of patho- 
genicity. 

Here w^e present a summary classification and a 
set of features for K massiliensis sp. nov. strain 
JG30T together v\^ith the description of the com- 
plete sequencing and annotation of its genome. 
These characteristics support the circumscription 
of the species K massiliensis. 

Classification and features 

A stool sample w^as collected from a healthy 16- 
year-old male Senegalese volunteer patient living 
in Dielmo [a rural village in the Guinean-Sudanian 
zone in Senegal), w^ho w^as included in a research 
protocol. The patient gave an informed and signed 
consent, and the agreement of the National Ethics 
Committee of Senegal and the local ethics commit- 
tee of the IFR48 [Marseille, France) were obtained 
under agreement 09-022. The fecal specimen w^as 
preserved at -80°C after collection and sent to 
Marseille. Strain JC30 [Table 1) w^as isolated in 
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January 2011 by aerobic cultivation on 5% sheep 
blood-enriched Columbia agar (BioMerieux). This 
strain exhibited a 96.9% nucleotide sequence sim- 
ilarity with K. gibsonii, the phylogenetically closest 
validated Kurthia species [Figure 1). This value 
was lower than the 97% 16S rRNA gene sequence 



threshold to delineate a new species without car- 
rying out DNA-DNA hybridization recommended 
by the report of the ad hoc committee on reconcil- 
iation of approaches to bacterial systematics [2]. 
Stackebrandt and Ebers proposed to increase this 
value to 98.7% [21]. 



Table 1. Classification and general features of Kurthia massiliensis strain JC30^ according to the 
MIGS recommendations [9] 



MIGS IP Property 



Term 



Evidence code^ 



Domain Bacteria 
Phylum Firmicutes 



IAS [10] 
IAS [11,12] 







Class Bacilli 


TAS 




Current classification 


Order Bacillales 


TAS 






Family Planococcaceae 


TAS 






Genus Kurthia 


TAS 






Species Kurthia massiliensis 


IDA 






Type strain JC30^ 


IDA 




Gram stain 


Positive 


IDA 




Cell shape 


Coccobacilli 


IDA 




Motility 


Motile by peritrichous flagella 


IDA 




Sporulation 


Nonsporulating 


IDA 




Temperature range 


Mesophile 


IDA 




Optimum temperature 


37°C 


IDA 


MIGS-6.3 


Salinity 


Growth in BHI medium + 3% NaCI 


IDA 


MIGS-22 


Oxygen requirement 


Aerobic 


IDA 




Carbon source 


Unknown 


NAS 




Energy source 


Unknown 


NAS 


MIGS-6 


Habitat 


Human gut 


IDA 


MIGS-15 


Biotic relationship 


Free living 


IDA 




Pathogenicity 


Unknown 


NAS 




Biosafety level 


2 




MIGS-14 


Isolation 


Human feces 




MIGS-4 


Geographic location 


Senegal 


IDA 


MIGS-5 


Sample collection time 


September 2010 


IDA 


MIGS-4.1 


Latitude 


13.7167 


IDA 


MIGS-4.1 


Longitude 


16.4167 


IDA 


MIGS-4.3 


Depth 


Surface 


IDA 


MIGS-4.4 


Altitude 


51 m above sea level 


IDA 



Evidence codes - IDA: Inferred from Direct Assay; TAS: Traceable Author Statement (i.e., a direct 
report exists in the literature); NAS: Non-traceable Author Statement (i.e., not directly observed for 
the living, isolated sample, but based on a generally accepted property for the species, or anecdo- 
tal evidence). These evidence codes are from the Gene Ontology project [19]. If the evidence is 
IDA, then the property was directly observed for a live isolate by one of the authors or an expert 
mentioned in the acknowledgements. 
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Kurthia zopfii DSM 20580' 

Kurthia sibsonii DSM 20636' 
Kurthia sihirica DSM 4747 ' 



Kurthia massiliensis strain JC30* 

Solibacillus silvestris strain HR3-23 



0.005 



Figure 1. Phylogenetic tree highlighting the position of Kurthia massiliensis strain JC30^ relative to other type strains 
within the genus Kurthia. Sequences were aligned using CLUSTALX, and phylogenetic inferences obtained using 
the neighbor-joining method within the MEGA 5 package [20]. Numbers at the nodes are percentages of bootstrap 
values obtained by repeating the analysis 1,000 times to generate a majority consensus tree. Solibacillus silvestris 
was used as outgroup. The scale bar represents 0.005 nucleotide change per nucleotide position. 



Surface colonies w^ere observed on sheep blood 
agar [bioMerieux) after 24 h aerobic incubation at 
37°C. The colonies of strain ]C30^ were circular, 
greyish/yellow^ish, shiny, curved and smooth, 2-5 
mm in diameter. Gram staining show^ed Gram- 
positive coccobacilli [Figure 2). 

Different grow^th temperatures [25, 30, 37, 45, 50 
and 55°C) w^ere tested. Grow^th occurred between 
25°C and 55°C, and optimal grow^th w^as observed 
betw^een 25°C and 50°C. Grow^th of the strain w^as 
tested under aerobic atmosphere, in the presence 
of 5% CO2, and under anaerobic and 
microaerophilic atmospheres, w^hich w^ere created 
using GENbag anaer and GENbag microaer 
[bioMerieux), respectively. The strains were aero- 
bic but also grew under microaerophilic condi- 
tions and in the presence of 5% CO2. Growth does 
not occur under anaerobic conditions. NaCl toler- 
ance of strain ]C30^ was determined on 
Difco™Brain Heart Infusion Agar plates [Becton 
Dickinson). The powder was supplemented with 
NaCl [Euromedex) to obtain the tested concentra- 
tions [0.5, 1, 2, 3, 5 10, 15%, w/v). Growth oc- 
curred between 0.5-5% NaCl but the optimum 
growth was between 0.5-3% NaCl. Growth in the 
range of pH 5.0-10.0 was tested using BEL™ Brain 
Heart Infusion [Becton Dickinson). pH tolerance 
revealed that growth could occur over a range of 
pH 6.0 - 9.0 with optimal growth between pH 7.0 - 
9.0. 

The size and ultrastructure of cells were deter- 
mined by negative staining transmission electron 
microscopy. The rods were 0.9-2.4 |im long and 
0.6-1.8 |j.m wide [Figure 3). Peritrichous flagella 
were observed. Capsule presence was determined 
by India ink stain and after bacteria embedding in 
Epon 812 resin and observation by transmission 



electron microscopy [Figures 4 and 5). Strain 
]C30^ exhibited catalase activity but no oxidase 
activity. Api ZYM, Api 20NE [BioMerieux) were 
used to study biochemical characters [Table 2]. 

Analysis of respiratory quinones by HPLC was 
carried out by the Identification Service and Dr 
Brian Tindall, DSMZ, Braunschweig, Germany. 
Respiratory lipoquinones were extracted from 
100 mg of freeze dried cell material as described 
by Tindall [22,23]. Respiratory lipoquinones were 
separated into their different classes 
[menaquinones and ubiquinones) by thin layer 
chromatography on silica gel, using hexane:ter- 
butylmethyl ether [9:1 v/v) as solvent. UV absorb- 
ing bands corresponding to menaquinones or 
ubiquinones were removed from the plate and 
further analyzed by HPLC with detection at 269 
nm. The only respiratory quinone for strain ]C30^ 
was MK-7 [100%). Preparation and determination 
of cellular fatty acids were carried out by follow- 
ing the procedures given for the Sherlock Microbi- 
al identification System [MIDI). The major fatty 
acids were Ci5:o iso 68.04% and Ci5:o anteiso 
16.92%. Polar lipids were extracted from 100 mg 
of freeze dried cell material using a chloro- 
form:methanol:0.3% aqueous NaCl mixture 
1:2:0.8 [v/v/v) [modified after [24]). The extrac- 
tion solvent was stirred overnight and the cell de- 
bris pelleted by centrifugation. Polar lipids were 
recovered into the chloroform phase by adjusting 
the chloroform:methanol:0.3% aqueous NaCl mix- 
ture to a ratio of 1:1:0.9 [v/v/v). Polar lipids were 
separated as previously described [25]. The polar 
lipids present were diphosphatidylglycerol, 
phosphatidylglycerol, phosphatidylethanolamine, 
phospholipid 1. The peptidoglycan of strain ]C30^ 
was isolated as described by Schleifer [26]. 
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Analysis was carried out as previously described 
[26,27] with the modification that TLC on cellu- 
lose was used rather than paper chromatography. 
Quantitative analysis of amino acids was per- 
formed following derivatization by gas chroma- 
tography and gas chromatography / mass spec- 
trometry (320-MS Quadrupole GC/MS, Varian) 
[28]. K. massiliensis showed the peptidoglycan 
type A4aL-Lys<^D-Glu (type A11.33 according to 
reference [36] ). 




Figure 2. Gram staining of K. massiliensis strain JC30^ 




Figure 3. Transmission electron microscopy of K. massiliensis 
strain JC30v using a Morgani 268D (Philips) at an operating 
voltage of GOkV.The scale bar represents 1 ^im. 



K. massiliensis was susceptible to penicillin 
amoxicillin, amoxicillin + clavulanic acid, 
imipenem, gentamycin, erythromycin, doxycy- 
cline, rifampicin, vancomycin, and nitrofurantoin. 
The organism was resistant to ceftriaxone, ciprof- 
loxacin, sulfamethoxazole trimethoprim and met- 
ronidazole. 
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Matrix-assisted laser-desorption/ionization time- 
of-flight [MALDI-TOF) MS protein analysis was car- 
ried out. Briefly, a pipette tip was used to pick one 
isolated bacterial colony from a culture agar plate, 
and to spread it as a thin film on a MALDI-TOF tar- 
get plate [Bruker Daltonics). Twelve distinct depos- 
its were made for strain ]C30^ from twelve isolated 
colonies and the manipulation was repeated anoth- 
er day. After air-drying, 1.5 |il matrix solution [sat- 
urated solution of a-cyanohydroxycinnaminic acid 
in 50% aqueous acetonitrile containing 2.5% 
trifluoroacetic acid) per spot was applied. MALDI- 
TOF MS was conducted using the Microflex LT 
spectrometer [Bruker Daltonics). All spectra were 
recorded in linear, positive ion mode. The accelera- 
tion voltage was 20 kV. Spectra were collected as a 
sum of 240 shots across a spot. Preprocessing and 
identification steps were performed using the 
manufacturer's parameters. The ]C30^ spectra 
were imported into the MALDI BioTyper software 
[version 3.0, Bruker) and analyzed by standard pat- 
tern matching [with default parameter settings) 
against the main spectra of 4,108 bacteria including 
those from K. gibsonii, K. sibirica and K. zopfii, used 
as reference data, in the BioTyper database. A score 
enabled the identification, or not, from the tested 
species: a score > 2.3 with a validly published spe- 
cies enabled the identification at the species level, a 
score > 1.7 but < 2 enabled the identification at the 
genus level; and a score < 1.7 did not enable any 
identification. For strain ]C30^, none of the ob- 
tained scores was > 1, thus suggesting that our iso- 
late was not a member of a known species. We in- 
cremented our database with the spectrum from 
strain ]C30^ [Figure 6). The spectrum was made 
available online in our free-access URMS database 
[29]. 




Figure 4. India ink capsule stain of K. massiliensis 
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Figure 5. Capsule characterization of K. massiliensis after the bacteria were 
embedded in Epon 812 resin and observed by transmission electron micros- 
copy. 



Table 2. Differential phenotypic characteristics between strain JC30^ 


and related 


species 


Characteristic 


1 


2 


3 


4 


gelatin hydrolysis 


+ 








N-acetyl-glucosamine assimilation 






+ 




D-maltose assimilation 


+ 








potassium gluconate assimilation 


+ 








trisodium citrate assimilation 


+ 








alkaline phosphatase 




+ 


w 


+ 


esterase (C4) 


+ 


+ 


w 


w 


esterase lipase (C8) 


+ 


+ 


w 


w 


valine arylaminidase 


w 




+ 




cystine arylaminidase 


+ 


+ 






Gchemotrypsin 


w 




+ 




naphtol-AS-BI-phosphohydrolase 






+ 


+ 


Gglucosidase 


+ 









Strains: 1, K. massiliensis sp. nov. JC30^; 2, K. gibsonii DSM 20636^; 3, K. zopfii DSM 
20580"^; 4, K. sibirica DSM 4747\ 

+: positive result, -: negative result, w: weak positive result 
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Figure 6. Reference mass spectrum from K. massiliensis strain JC30^. Spectra from 24 individual colonies 
were compared and a reference spectrum was generated. 



Genome sequencing information 

Genome project history 

The organism was selected for sequencing on the 
basis of its phylogenetic position and 16S rRNA 
similarity to other members of the genus Kurthia, 
and is part of a "culturomics" study of the human 
digestive flora aiming at isolating all bacterial spe- 
cies w^ithin human feces. It w^as the first genome of 
a Kurthia species A summary of the project infor- 
mation is show^n in Table 3. The EMBL accession 
number is CAEUOIOOOOOO and consists of 98 
contigs (>200 bp) and 18 scaffold (> 2,424 bp). 
Table 3 show^s the project information and its as- 
sociation w^ith MIGS version 2.0 identifiers. 

Growth conditions and DNA isolation 

K massiliensis sp. nov. strain ]C30^, CSUR P141T 
DSM 24639'^; was grown aerobically on 5% sheep 
blood-enriched Columbia agar at 37°C. Three petri 
dishes were spread and resuspended in 3x100 |j.l 
of G2 buffer. A first mechanical lysis was per- 
formed by glass powder on the Fastprep-24 de- 
vice [Sample Preparation system) from MP 
Biomedicals, USA using 2x20 second cycles. DNA 
was then treated with lysozyme [4.17g/L, 30 
minutes at 37°C) and extracted through the 
BioRobot EZ 1 Advanced XL (Qiagen). The DNA 



was then concentrated and purified on a Qiamp kit 
[Qiagen). The yield and the concentration were 
measured by the Quant-it Picogreen kit [Invitro- 
gen) on the Genios Tecan fluorometer at 63.1/|j.l. 

Genome sequencing and assembly 

Shotgun and 3-kb paired-end sequencing strategies 
were used. The shotgun library was constructed 
with 500 ng of DNA with the GS Rapid library Prep 
kit [Roche). For paired-end sequencing, 5 |ig of DNA 
was mechanically fragmented on a Hydroshear de- 
vice [Digilab) with an enrichment size at 3-4 kb. The 
DNA fragmentation was visualized using the 2100 
BioAnalyzer [Agilent) on a DNA labchip 7500 with 
an optimal size of 3.619 kb. The library was con- 
structed according to the 454 GS FLX Titanium 
paired-end protocol. Circularization and nebuliza- 
tion were performed and generated a pattern with 
an optimal size of 472 bp. After PGR amplification 
through 15 cycles followed by double size selection, 
the single stranded paired-end library was then 
quantified using the Genios fluorometer [Tecan) at 
430 pg/|iL. The library concentration equivalence 
was calculated as 1.69E+09 molecules/|iL. The li- 
brary was stored at -20°C until further use. 
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Table 3. Project information 



MIGS ID 


Property 


Term 


MIGS-31 


Finishing quality 


High-quality draft 


MIGS-28 


Libraries used 


One paired end 3-kb library and one Shotgun library 


MIGS-29 


Sequencing platforms 


454 GS FLX Titanium 


MIGS-31. 2 


Fold coverage 


22x 


MIGS-30 


Assemblers 


Newbler version 2.5.3 


MIGS-32 


Gene calling method 


Prodigal 




EMBL ID 


CAEUOl 000000 




EMBL Date of Release 


February 12, 2012 




Project relevance 


Study of the human gut microbiome 



The shotgun and paired-end libraries were cfonally- 
amphfied with 3 cpb and Icpb in 3 and 4 emPCR 
reactions respectively on the GS Titanium SV emPCR 
Kit [Lib-L) v2 [Roche). The yields of the emPCR were 
18.65 and 14.31% respectively. Approximately 
340;000 beads for the shotgun sequencing and 
790,000 beads for the 3kb paired end sequencing 
were loaded onto the GS Titanium PicoTiterPlate 
PTP Kit 70x75 and sequenced with the GS FLX Tita- 
nium Sequencing Kit XLR70 [Roche). The run was 
performed overnight and then analyzed on the clus- 
ter through the gsRunBrowser and Newbler assem- 
bler [Roche). A total of 294,263 passed filter wells 
were obtained and generated 81.3 Mb with a length 
average of 301 bp. The passed filter sequences were 
assembled using Newbler with 90% identity and 40 
bp as overlap. The final assembly identified 18 scaf- 
folds and 72 large contigs [> 1,500 bp). 

Genome annotation 

Coding sequences [CDSs) were predicted using 
PRODIGAL with default parameters [30]. The func- 
tional annotation of protein sequences was per- 



formed against the non-redundant GenBank data- 
base using BLASTP. Functional categories of these 
proteins were searched against the Clusters of 
Orthologous Groups [COG) database using 
COGNITOR [31]. The prediction of RNAs genes, i.e., 
rRNAs, tRNAs and other RNAs was carried out using 
RNAmmer [32] and ARAGORN [33] algorithms. The 
transmembrane segments and peptide signals were 
identified using TMHMM [34] and SignalP tools [35]. 

Genome properties 

The genome is 3,199,090 bp long with a 39.26% 
GC content [Table 4, Figure 7). Of the 3,326 pre- 
dicted genes, 3,240 were protein-coding genes, 
and 86 were RNAs. A total of 2,425 genes [74.8%) 
were assigned a putative function. The remaining 
genes were annotated as either hypothetical pro- 
teins or proteins of unknown functions. The dis- 
tribution of genes into COGs functional categories 
is presented in Table 5. The properties and the 
statistics of the genome are summarized in Tables 
4 and 5. 



Table 4. Nucleotide content and gene count levels of the genome 



Attribute 


Value 


% of total 


Genome size (bp) 


3,199,090 




DNA coding region (bp) 


2,794,828 


87.4 


DNA G+C content (bp) 


1,255,894 


39.26 


Total genes 


3,326 


100 


RNA genes 


86 


2.6 


Protein-coding genes 


3,240 


97.4 


Genes with function prediction 


2,425 


74.8 


Genes assigned to COGs 


2,500 


77.16 


Genes with peptide signals 


269 


8.3 


Genes with transmembrane helices 


473 


14.6 



a) The total is based on either the size of the genome in base pairs or the 
total number of protein coding genes in the annotated genome. 



http://standardsingenomics.org 



227 



Kurthia massiliensis sp. nov. 




Figure 7. Graphical circular map of Kurthia massiliensis genome. From outside to the center: Genes on both 
strands, genes on foward strand, genes on reverse strand and genes colored by COG categories. 
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Table 5. Number of genes associated with the 25 general COG functional categories 
Code Value %age Description 



J 


161 


4.97 


Translation 


A 


0 


0 


RNA processing and modification 


K 


218 


6.73 


Transcription 


L 


184 


5.68 


Replication, recombination and repair 


B 


1 


0.03 


Chromatin structure and dynamics 


D 


34 


1.05 


Cell cycle control, mitosis and meiosis 


Y 


0 


0 


Nuclear structure 


V 


47 


1.45 


Defense mechanisms 


T 


171 


5.28 


Signal transduction mechanisms 


M 


118 


3.64 


Cell wall/membrane biogenesis 


N 


82 


2.53 


Cell motility 


Z 


0 


0 


Cytoskeleton 


w 


0 


0 


Extracellular structures 


u 


43 


1.33 


Intracellular trafficking and secretion 


o 


92 


2.84 


Posttranslational modification, protein turnover, chaperones 


c 


139 


4.29 


Energy production and conversion 


G 


134 


4.14 


Carbohydrate transport and metabolism 


E 


267 


8.24 


Amino acid transport and metabolism 


F 


74 


2.28 


Nucleotide transport and metabolism 


H 


134 


4.14 


Coenzyme transport and metabolism 


1 


110 


3.40 


Lipid transport and metabolism 


P 


204 


6.30 


Inorganic ion transport and metabolism 


Q 


68 


2.10 


Secondary metabolites biosynthesis, transport and catabolism 


R 


402 


12.41 


General function prediction only 


S 


237 


7.31 


Function unknown 


X 


740 


22.84 


Not in COGs 



The total is based on the total number of protein coding genes in the annotated genome. 
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Comparison with other Kurthia genomes 

To date, no genome of other strains or species be- 
longing to the genus Kurthia were sequenced. 

Conclusion 

On the basis of phenotypic, phylogenetic and ge- 
nomic analyses, we formally propose the creation 
of Kurthia massiUensis sp. nov., which contains the 
strain JCSO'^. This bacterium was found in Senegal. 

Description of Kurthia massiliensis sp. nov. 

Kurthia massiUensis [mas.si.li.en'sis. L. masc. adj. 
massiliensis of Massilia, the old Roman name for 
Marseille, where the type strain was isolated). Iso- 
lated from stool of a healthy Senegalese patient. K 
massiUensis are aerobic Gram-positive coccobacilli. 
On sheep blood agar after 24 h aerobic incubation 
at 37°C, colonies of strain JCSO'^ are circular, grey- 
ish/yellowish, shiny, curved and smooth, 2-5 mm in 
diameter. Cells are motile by peritrichous flagella 
and encapsulated. Catalase activity is positive but 



oxidase activity is negative. Gelatine hydrolysis, 
maltose assimilation, potassium gluconate assimi- 
lation, malic acid assimilation, trisodium citrate 
assimilation are present. Esterase [C4), esterase 
lipase [C8), cystine arylaminidase, a-gluconidase 
activities are observed. Valine arylaminidase and 
alpha-chemotrypsin activities are weakly positive. 
The major fatty acids are iso Ci5:o 68.04% and 
anteiso C15:0 16.92%. Polar lipids found are 
diphosphatidylglycerol, phosphatidylglycerol, 
phosphatidylethanolamine, and phospholipid 1. 
The peptidoglycan type is A4aL-Lys^D-Glu (type 
All. 33 according to [36]). Cells are susceptible to 
penicillin G, amoxicillin, amoxicillin + clavulanic 
acid, imipenem, gentamycin, erythromycin, 
doxycycline, rifampicin, vancomycin and 
nitrofurantoin. The G+C content of the genome is 
39.26%. The type strain is JC30T (= CSUR P141T = 
DSM 24639T). 
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